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Periodic Phase Synchronization in coupled chaotic oscillators 



o 
o 

> 
O 

0\ 



Q 

U 
d 



> 

in 
o 
in 
o 
m 
O 



> 
X 



Won-Ho Kye 1 ^] Dae-Sic Lee 1,2 , Sunghwan Rim 1 , Chil-Min Kim 1 !]] and Young- Jai Park 2 
1 National Creative Research Initiative Center for Controlling Optical Chaos, 
Pai-Chai University, Daejeon 302-735, Korea and 
2 Department of Physics, Sogang University, Seoul 121-742, Korea 

We investigate the characteristics of temporal phase locking states observed in the route to phase 
synchronization. It is found that before phase synchronization there is a periodic phase synchro- 
nization state characterized by periodic appearance of temporal phase locking state and that the 
state leads to local negativeness in one of the vanishing Lyapunov exponents. By taking a statistical 
measure, we present the evidences of the phenomenon in unidirectionally and mutually coupled 
chaotic oscillators, respectively. And it is qualitatively discussed that the phenomenon is described 
by a nonuniform oscillator model in the presence of noise. 
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Over the past decade, synchronization in chaotic oscil- 
lators P, has received much attention because of its 
fundamental importance in nonlinear dynamics and po- 
tential applications to laser dynamics [3(, electronic cir- 
cuits Q, chemical and biological systems and secure 
communications 6]. Synchronization in chaotic oscilla- 
tors is characterized by the loss of exponential instabil- 
ity in the transverse direction through interaction. In 
coupled chaotic oscillators, it is known, various types of 
synchronization are possible to observe, among which are 
complete synchronization (CS) P, , phase synchroniza- 
tion (PS) |7j,|8|], lag synchronization (LS) [9( and gener- 
alized synchronization (GS) jlCj . 

One of the noteworthy synchronization phenomena in 
this regard is PS which is defined by the phase lock- 
ing between nonidentical chaotic oscillators whose ampli- 
tudes remain chaotic and uncorrelated with each other: 
|#i — 6>2 1 < const. Since the first observation of PS in 
mutually coupled chaotic oscillators Q, there have been 
extensive studies in theory Q and experiments jj. The 
most interesting recent development in this regard is 
the report that the interdependence between physiologi- 
cal systems is represented by PS and temporary phase- 
locking (TPL) states, e.g., (a) human heart beat and res- 
piration [ill , (b) a certain brain area and the tremor ac- 
tivity |l2j|. etc 0,^3- Application of the concept of PS 
in these areas sheds light on the analysis of nonstation- 
ary bivariate data coming from biological systems which 
was thought to be impossible in the conventional statis- 
tical approach. And this calls new attention to the PS 
phenomenon. 

Accordingly, it is quite important to elucidate a de- 
tailed transition route to PS in consideration of the recent 
observation of a TPL state in biological systems. What is 
known at present is that TPLQ transits to PS and then 
transits to LS as the coupling strength increases. On the 
other hand, it is noticeable that the phenomenon from 



nonsynchronization to PS have hardly been studied, in 
contrast to the wide observations of the TPL states in 
the biological systems. 

The chief goal of this Letter is to study the char- 
acteristics of TPL states observed in the regime from 
nonsynchronization to PS in coupled chaotic oscillators. 
We report that there exists a special locking regime in 
which a TPL state shows maximal periodicity, which phe- 
nomenon we would call periodic phase synchronization 
(PPS). We show this PPS state leads to local negative- 
ness in one of the vanishing Lyapunov exponents, taking 
the measure by which we can identify the maximal peri- 
odicity in a TPL state. We present a qualitative expla- 
nation of the phenomenon with a nonuniform oscillator 
model in the presence of noise. 

We consider here the unidirectionally coupled non- 
identical Rossler oscillators for first example: 



x\ = -u>iy\ - zi, 

yi = uiixi + 0.15j/i, 

zi = 0.2 + zi(xi - 10.0), 

±2 = -UJ2V2 ~ Zl, 

y 2 = uj 2 x 2 + 0.165y 2 + e(yi 

Z2 = 0.2 + z 2 {x 2 - 10.0), 
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where the subscripts imply the oscillators 1 and 2, re- 
spectively, CJ1.2 (= 1.0 ± 0.015) is the overall frequency 
of each oscillator, and e is the coupling strength. It is 
known that PS appears in the regime e > e c and that 
2ir phase jumps arise when e < e c . Lyapunov exponents 
play an essential role in the investigation of the transi- 
tion phenomenon with coupled chaotic oscillators and as 
generally understood that PS transition is closely related 
to the transition to the negative value in one of the van- 
ishing Lyapunov exponents 0- 

Figure 1 shows two largest conditional Lyapunov expo- 
nents from Eq. (1) according to the coupling strength e. 
One can see the dip characterized by the local negative- 
ness in the vanishing Lyapunov exponent. The reference 
points A and C indicate the borders of the dip and B its 
center. PS transition occurs in the right of C (e = 0.085) 
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FIG. 1: Two largest conditional Lyapunov exponents and the 
dip in coupled Rossler oscillators when u>o — 1.0: A (e = 
0.058), B(e = 0.072), and C (e = 0.083) are the reference 
points to explain different characteristics of phase dynamics. 



where the phase difference of coupled chaotic oscillators 
is bounded within a constant. The temporal behaviors 
at the three reference points are presented in Fig. 2. The 
phase jumps at B looks quite regular, compared to those 
of A and C. Though this observation may seem rather 
intuitive, we shall see that it is a valid one and that all 
the phenomenon is deeply related to the dip, the local 
negativeness in the vanishing Lyapunov exponent of Fig. 
1. 

A vanishing Lyapunov exponent corresponds to a 
phase variable of an oscillator and it exhibits the neutral- 
ity of an oscillator in the phase direction. Accordingly, 
the local negativeness of an exponent indicates this neu- 
trality is locally broken 0. It is important to define an 
appropriate phase variable in order to study the TPL 
state more thoroughly. In this regard, several methods 
have been proposed methods of using linear interpola- 
tion at a Poincare section pfl, phase space projection 
00 , tracing of the center of rotation in phase space [l5| . 
Hilbert transformation ^| , or wavelet transformation 
[l4|. Among these we take the method of phase space 
projection onto the x\ — y\ and x-i — yi planes with the 
geometrical relation 0± t 2 — ajctan.($/i a/a^i 2)) and obtain 
phase difference <j> = 6\ — Q^. 

The system of coupled oscillators is said to be in a TPL 
state (or laminar state) when (<f>) < A c where (...) is the 
running average over appropriate short time scale and 
A c is the cutoff value to define a TPL state. The locking 
length of the TPL state, r, is defined by time interval 
between two adjacent peaks of ((f)) (see Fig. 2). In order 
to study the characteristics of the locking length r, we 
introduce a measure: 
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FIG. 2: Temporal behaviors of the phase difference: (a) <f> at 
each reference point, r is the locking length of a TPL state, 
(b) {(f>) at reference point C (running averaged over 10 seconds 
to remove irrelevant fast fluctuations). The cutoff value for 
TPL a state is A c = 1.0 x 10" 4 . 
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FIG. 3: Measure P as a function of e. Big dots are the ref- 
erence points of Fig.l and each point is the average value of 
2500 TPL states. 



which is the ratio between the average value of time 
lengths of TPL states and their standard deviation. In 
terminology of stochastic resonance, it can be interpreted 
as noise-to-signal ratio 0, . The measure would be 
minimized where the periodicity is maximized in TPL 
states. 

Measure P as a function of e is presented in Fig. 3. The 
big dots indicate the reference points of Fig. 1 . We see 
that the value of P begins to drop rapidly from reference 
point A which corresponds to the left border of the dip 
in Fig. 1. And the value of P is minimized in a broad 
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region around B, the center of the dip. The value rapidly 
increases after passing reference point C, the right border 
of the dip. What is interesting here is that in the region 
from e = 0.062 to 0.078 the periodicity is maximized and 
corresponds to the central part of the dip. Eventually, 
the coupled chaotic oscillators develop to PS near e = 
0.085. The result presented in Fig. 3 leads us to argue 
that the dip in a vanishing Lyapunov exponent shows the 
maximal periodicity of TPL states around the minimum 
of P. We call the TPL state inside the dip a PPS state in 
the sense that the TPL state appears rather periodically 
than outside the dip. 

To validate the argument, we explain the phenomenon 
in simplified dynamics. From Eq. (1), we obtain the 
equation of motion in terms of phase difference: 

^ = ALu + A{9 1 ,e 2l e)smcj ) + ^e u 6 2 ,e), (3) 

where, 

A(6 1 ,0 2 ,e) = (£ + 0.15)008(0! + 2 )-^ A, 

l ti 2 

z n 2 Hi K 2 

+ (e + 0.015) cos((9 2 )sin((92). 

Here Auj = uj\ — uj 2 and R\ 2 = \J x\ 2 + y\ 2 - And from 

Eq. (3) we obtain the simplified equation to describe the 
phase dynamics: d<fi/dt — Auj + (A) sin((/>) + £, where 
(A) is the time average of A{6\, 8 2 ,e)- This is a nonuni- 
form oscillator in the presence of noise where £ plays 
a role of effective noise [l(| and the value of (A) con- 
trols the width of bottleneck (i.e, nonuniformity of the 
flow). If the bottleneck is wide enough, (i.e., faraway 
from the saddle- node bifurcation point: Auj 3> —(A)), 
the effective noise hardly contributes to the phase dy- 
namics of the system. So the passage time is wholly 
governed by the width of the bottleneck as follows: 
(r) ~ I/^/Auj 2 - (A) 2 ~ 1/v/Aw 2 -e 2 /4, which is a 
slowly increasing function of e. In this region while the 
standard deviation of TPL states is nearly constant (be- 
cause the widely opened bottlenecks periodically appears 
and those lead to small standard deviation), the average 
value of locking length of TPL states is relatively short 
and the ratio between them is still large. Accordingly, 
the value of P(e) slowly decreases in the regime before 
reference point B in Fig. 3. 

On the contrary as the bottleneck becomes narrower 
(i.e., near the saddle-node bifurcation point: Auj > 
— (A)) the effective noise begins to perturb the process 
of bottleneck passage and regular TPL states develop 
into intermittent ones (see C in Fig. 2 (a)) |a, LiJ- It 
makes the standard deviation increase very rapidly and 
this trend overpowers that of the average value of locking 
lengths of the TPL states. For that reason, the value of 
P(e) rapidly increases passing the PPS regime in Fig. 3. 
Thus we understand that the competition between width 
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FIG. 4: Lyapunov exponents and measure P in mutually 
coupled Rossler oscillator: ±1,2 = — 0/1,23/1,2 — 21,2 + e(x2,i — 

3:1,2), Vl,2 = Ul,2Xl,2 + 0.151/1,2,21,2 = 0.2 + Zipfalfi - 10.0), 

where 0/1,2 = 1.0 ± 0.015. (a) One can see the dip in one of 
the vanishing Lyapunov exponents and the reference point D 
indicates the center of it. (b) The big dot shows the minimum 
point and it coincides with the point of the center of the dip. 
Each point is the average value of 12000 TPL states. 



of bottleneck and amplitude of effective noise produces 
the crossover at the minimum point of P(e) which shows 
the maximal periodicity of TPL states. 

Rosenblum et al. firstly observed the dip in mutu- 
ally coupled chaotic oscillators 0. However the origin 
and the dynamical characteristics of the dip have been 
left unclarified. We argue that the dip observed in mu- 
tually coupled chaotic oscillators has the same origin as 
observed above in unidirectionally coupled systems. Fig- 
ure 4 shows the four largest Lyapunov exponents and the 
value of P according to the coupling strength e. We can 
see that the minimum region of P from e = 0.022 to 
0.027 coincides with the central part of the dip (D in the 
figure) even though the valley in Fig. 4 (b) is not deeper 
than that of the unidirectionally coupled systems. Thus 
it is reasonably summed up that a PPS state exists just 
before the transition to PS and that the dip observed by 
Rosenblum et al. is the very evidence of a PPS state in 
mutually coupled chaotic oscillators. 

Common apprehension is that near the border of syn- 
chronization the phase difference in coupled regular os- 
cillators is periodic whereas in coupled chaotic os- 
cillators it is irregular Q. On the contrary, we report 
that the special locking regime exhibiting the maximal 
periodicity of a TPL state also exists in the case of cou- 
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pled chaotic oscillators. In general, the phase difference 
of coupled chaotic oscillators is described by the one- 
dimensional Langevin equation: <fi — F((j)) + £ where £ is 
the effective noise with finite amplitude. The investiga- 
tion with regard to PS transition is the study of scaling 
of the laminar length around the virtual fixed point <p* 
where F{<j>*) = [H HI and PS transition is estab- 
lished when \ f£ F((f>)d(f>\ > max|£|. Consequently, the 
crossover region, from which the value of P grows ex- 
ponentially (as shown in Fig. 3-4), exists because inter- 
mittent series of TPL states with longer locking length r 
appears as PS transition is nearer. Eventually it leads to 
an exponential growth of the standard deviation of the 
locking length. Thus we argue that PPS is the generic 
phenomenon mostly observed in coupled chaotic oscilla- 
tors prior to PS transition. 

In conclusion, analyzing the dynamic behaviors in cou- 
pled chaotic oscillators with slight parameter mismatch 
we have completed the whole transition route to PS. 
We find that there exists a special locking regime called 
PPS in which a TPL state shows maximal periodicity 
and that the periodicity leads to local negativeness in 
one of the vanishing Lyapunov exponents. We have also 



made a qualitative description of this phenomenon with 
the nonuniform oscillator model in the presence of noise. 
Investigating the characteristics of TPL states between 
nonsynchronization and PS, we have clarified the tran- 
sition route before PS. Since PPS appears in the inter- 
mediate regime between non-synchronization and PS, we 
expect that the concept of PPS can be used as a tool for 
analyzing weak interdependences, i.e. those not strong 
enough to develop to PS, between nonstationary bivari- 
ate data coming from biological systems, for instance. 
Moreover PPS could be a possible mechanism of the 
chaos regularization phenomenon 0, [2^ observed in 
neurobiological experiments. 

Note added. - Recently, we were informed by S. Boc- 
caletti that the phenomenon observed by us was con- 
firmed in CO2 laser systems, experimentally j^. 
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about experimental confirmation of PPS and thank A. 
Pikovsky, Y. -C. Lai, K. Josic, and M. Choi for helpful 
comments. This work is supported by Creative Research 
Initiatives of the Korean Ministry of Science and Tech- 
nology. 
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